Quantum fluctuations in thin superconducting wires of finite length 
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In one dimensional wires, fluctuations destroy superconducting long-range order and stiffness 
at finite temperatures; in an infinite wire, quasi-long range order and stiffness survive at zero 
temperature if the wire's dimensionless admittance fi is large, fi > 2. We analyze the disappearance 
of this superconductor-insulator quantum phase transition in a finite wire and its resurrection due 
to the wire's coupling to its environment characterized through the dimensionless conductance K. 
Integrating over phase slips, we determine the flow of couplings and establish the \i-K phase diagram. 
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In systems with reduced dimensions (films and wires) 
fluctuations and disorder strongly influence on the super- 
conducting transition temperature []J , eventually driving 
a superconductor-insulator (SI) transition which has at- 
tracted a lot of attention recently 0, EJ ■ In ID wires, 
the fluctuations of Cooper pairs appearing below the 
mean-field transition at T c o define a finite resistance via 
nucleation of thermally activated phase slips [f| and 
hence remove the finite temperature transition; super- 
conductivity possibly survives only at zero temperature. 
The focus then is on the quantum nucleation of phase 
slips; their proliferation may trigger a zero temperature 
SI quantum phase transition IE S EL While first 
attempts to observe quantum phase slips 7| are still de- 
bated due to the granular structure of the wires [llLJ l3 l . 
recent experiments on amorphous ultra-thin wires Li 
carry the signatures of a SI transition in a homogeneous 
system. In this letter, we analyze the bosonic Cooper 
pair fluctuations in realistic wires of finite length which 
are naturally coupled to their environment through their 
boundaries; we demonstrate how the SI transition is 
quenched in the finite system and reappears through its 
coupling to the environment. 

Previous studies of bosonic fluctuations in infinite 
wires Q have found a T — SI quantum phase tran- 
sition at the critical value /i = 2 of the dimensionless ad- 
mittance characterizing the wire's superconducting prop- 
erties. At T > the system exhibits a finite resistance; 
ignoring an additional dissipative channel due to excited 
quasi-particles, we call the non-superconducting state an 
insulating one. On the other hand, real experiments are 
carried out on wires of finite length, L ~ 0.5 — 2 /im 
typically The wire's coupling to the environment 
through appropriate boundary conditions imposes a dras- 
tic change in the phase slip dynamics and modifies the 
wire's low-energy physics. A generic description is ob- 
tained by embedding the wire in a voltage driven loop 
with impedances Z p (uj) and Z s (lj) placed in parallel and 
in series, see Fig. ^ Here, we are mainly interested in 
thermodynamic aspects (involving the static resistances 
R p and R s ); furthermore we concentrate on the current- 
driven limit with R s ^S> R P , then / = V/R s is fixed. 
The voltage driven case with R p <C R s and other mixed 
cases (R s of order R p ) are easily derived from the current- 
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FIG. 1: Left: Setup with a quantum wire embedded in a 
voltage (V) driven loop with parallel (Z p ) and serial (Z B ) 
impedances defining the environment. Right: Phase diagram 
with superconducting and insulating phases separated by a 
quantum phase transition at K = 1. The superconducting 
phase splits into weak and strong regimes separated by a 
crossover at /i ~ 2, the leftover of the SI transition in the 
infinite wire. The insets show sketches of the wire's I w -V v 
characteristic at T — 0. In the superconductor (K > I) the 
algebraic characteristic is dominated by the environment at 
small currents I < /dec- A highly conducting shunt with 
K 2> /i allows to probe the wire above the deconfinement 
current Idee', Ic denotes the critical current of the wire. The 
insulator exhibits a Coulomb gap behavior below the critical 
voltage V c - 



driven solution via Kirchhoff 's laws, see below. 

The coupling to the environment changes the T = 
phase diagram of the infinite system, see Fig. ^ The 
SI transition at /J, = 2 is turned into a crossover, while 
a new quantum phase transition appears at the criti- 
cal value K = Rq/R p = 1 of the shunt's conductance 
(Rq = nTi/2e 2 denotes the quantum resistance). A well 
conducting shunt with K > 1 relaxes the strain on the 
wire and produces a superconducting response with an 
algebraic I-V characteristic, while a low conductance 
K < 1 leads to the proliferation of phase slips and hence 
to an insulator. The characteristic response of the super- 
conducting wire can be probed with a high conductance 
environment K /i at large drive above the current /dec 
where confined phase-slip pairs are separated. 

In thin superconducting wires the plasma mode ac- 
quires a linear dispersion H3 | with sound velocity c s 
(this contrasts with the bulk case where the Coulomb 
interaction lifts the plasma mode to finite frequencies). 
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The low frequency action describing the bosonic fluctua- 
tions of the Cooper pairs then takes the form 0, 0, [^j 



S w 



2?r 



h/3 pL/2 

dT I dx 

J-L/2 



c s (d x (f)) 2 + — {d T <ff 

C.s 



(1) 



the first (second) term accounts for the kinetic (compres- 
sion) energy (we assume a high energy cutoff H/k limiting 
the validity of QJ). Within a mean-field description the 
dimensionless admittance \x and the sound velocity c s 
are related to the ID superfluid density p s — n s S/2m* 
and the Coulomb interaction between the Cooper pairs 
HG3, hfic s /2Tr = h 2 p s /2 and hfj,/2-n:c s = {h/2e) 2 C/2 
(here, S is the wire's cross section and C = e/ \n(d 2 / S) is 
the capacitance per unit length, with e the dielectric con- 
stant of the surrounding media placed a distance d away) . 
With Al = (47m s e 2 /rac 2 ) 1 / 2 the London penetration 
length we obtain /x » 30VCS/A L and c 2 Jc 2 w 0.1S/A|C. 
Going beyond the mean-field level, both fermionic and 
bosonic high-energy fluctuations play a critical role in 
the proper determination of the low-energy action |l5j : 
here, we assume such effects to be included in our choice 
of effective phenomenological parameters fi and c s . 

The coupling between the superconducting wire and 
the environment involves the boundary fields 4>±{t) = 
0±{x,t)\ x=l/2 where (/)±(x,t) = 0(x,t)±^(-x,t): fluc- 
tuations in the phase difference generate a voltage V w 
across the wire inducing currents in the parallel shunts, 
while fluctuations in <p + account for charge accumula- 
tion. Here, we concentrate on the current-driven setup 
(see Fig. ^) described by the low-energy action 
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with an ideal current source (R s >• R p ) driving the sys- 
tem with the current / and a parallel resistor with resis- 
tance R p accounting for the dissipation (K = Rq/R p is 
the dimensionless conductance of the shunt). Protecting 
the system from the measurement setup using appropri- 
ate filtering 0], we can neglect capacitive effects. The 
environment (J2J does not account for an intrinsic, e.g., 
quasi-particle induced, dissipation in the superconduct- 
ing wire; the latter contributes to the precise values of the 
parameters fx, c s , and A (the vortex fugacity, see below). 

The statistical mechanics of the system 'wire plus en- 
vironment' is determined by the partition function 



Z = J D[^\exp[-S(4>)/n} 



(3) 



with S = S w + S c . Its main contributions arise from the 
combination of Gaussian and topological fluctuations, 
so-called phase slips or instantons. Gaussian fluctua- 
tions destroy the superconducting long-range order as 
expressed by the logarithmically diverging phase corre- 
lator ([<f)(x, t) — cj)(x, t')] 2 ) ~ K^ 1 In |t — t'|; however, 



the surviving quasi-long range order is sufficient to allow 
for a finite phase stiffness, i.e., superconducting response. 
The latter is destroyed by the proliferation of phase slips, 
the process we are going to analyze in detail now. 

Quantum phase slips are vortex-like solutions in x, t- 
space with finite winding [y = ±1) around a core region 
of size x c < c s k and r c < k where the superconducting 
gap A(x,t) drops to zero; outside this core region they 
satisfy the differential equation [c 2 <9 2 + d^](p(T,x) = 
within the wire and respect the boundary conditions 



K\uj\4>- + fj,c s [d x 4i\- = nl/e 



M+ - 0. (4) 



The first equation describes current conservation as given 
by Kirchhoff's law, with the supercurrent I w = (e/ic s /7r) 
[<9a;<?!>]_ through the wire and the dissipative current 
V w I R p in the shunt adding up to the total external cur- 
rent / (here, [d x (j)]±(T) = d x (j)±(x, r)| x=i / 2 and V w = 
h\oj\(f>-/2e is the voltage over the wire). The second equa- 
tion accounts for charge neutrality; we ignore additional 
capacitive contributions here. Note that Kirchhoff's law 
maintains the form (0J when going over to the general 
situation with impedances both in series and in parallel, 
Z w + V w /R t = V/Rs with Rt 1 = R~ 1 + R- 1 the total 
conductance and the total current / = V/ R s is expressed 
through the driving voltage V. Hence once we know the 
solution V W (I) = Rqxi{I) to the current-driven problem 
(with Rp = R), the solution I W (V) = Xv(V/Rq) to the 
voltage-driven situation (with R s = R) follows from the 
relation / - Kxi(I) = Xv(I/K), with K = R Q /R. 

The partition function @ factorizes into Gaussian and 
topological parts, Z — Z G Z top , where the latter can be 
expanded in a series Z top = J2^ =0 {X n /n\) 2 Z n {G) with 
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describing a (neutral) gas of n vortex-anti- vortex pairs, 
each contributing with a dimensionless action G = 
Spni r /h (we assume 1 = 0): in addition, each vortex-pair 
is weighted with the fugacity A 2 accounting for the micro- 
scopic structure of the cores. Here, we assume that this 
factor is small, A ~ exp(—ARQL/ Rn£) <C 1 with Rn 
the normal resistance of a wire, £ the superconducting 
coherence length, and A a factor of order unity 0- 

The interaction G between vortex-anti-vortex pairs is 
the crucial quantity determining the system's behavior 
- here we only quote the relevant results, see Ref. |J3 
for details. For the infinite wire, the interaction between 
vortices 0k = ± arg[(x - xu) + i c s (t - Tk)], k = 1, 2, is 
logarithmic at all distances x — x± — X2, f = t\ — T2, 



G(x,t) = nln[{x 2 /c 2 s +t 2 )/k 2 



(5) 



A standard Kosterlitz-Thouless (KT) scaling analysis 
provides the RG equations = — 47r 2 /i 2 A 2 and 
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di\ = (2 — fx)X (I is the scaling parameter) and the sys- 
tem undergoes a Berezinskii-Kosterlitz-Thouless transi- 
tion 0| at /x = 2. For fi > 2, the number of free 
vortex-anti-vortex pairs is quenched and the wire is in 
the superconducting phase; for < 2 the vortex fugacity 
increases and free vortex-anti- vortex pairs proliferate — 
the system turns insulating. 

In a finite wire the solution has to respect the bound- 
ary condition ijljl. At low frequencies and vanishing 
drive these reduce to the Neumann boundary conditions, 
d x <fi N (±L/2,T) = 0, i.e., the currents induced by phase 
slips cannot leave the wire. The solution (j) K (x, r) de- 
rives from the 2L-periodic solution 4> 2 l- p using mirror 
vortices, N (x, r) = cf> 2L _ p (x,T) + (f> 2L _ p (L-x 1 T); in turn, 
the instanton solution with periodic boundary conditions 
4>2l- v {x,t) = 4> 2 L-p{x + 2L,t) derives from the solution 
in the infinite wire via conformal transformation. Fig. 
12 illustrates the geometry of the phase slips; most re- 
markably, at large distances the mutual screening of the 
defect-pair is replaced by the screening via image charges. 
Hence the wire's contribution G w to the action saturates 
at inter-vortex separations r = \t% — t± \ > L/nc s , result- 
ing in asymptotically free vortices. Furthermore, their 
action can be minimized by moving them to the bound- 
ary such that G w = 0. 
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FIG. 2: Phase slip solutions with periodic (solid lines denote 
contours 4> P — const., left) and Neumann (0 N , right) bound- 
ary conditions at intervortex distance L/c 3 < f. Left: for 
K 2> /i, defects are screened mutually, the 2-zr phase drop 
appears along the i-axis and drives a large current (dashed) 
through the highly conducting shunt; the resulting string con- 
fines defect pairs. Right: with /i K, defects are screened 
individually by their mirror images (grey underlaid), the 2n 
phase drop appears along the r-axis and sets up displacement 
currents in the wire which cannot escape through shunt, hence 
charge accumulates at the boundary resulting in a large volt- 
age over the shunt; the defects are asymptotically free. 

In this situation the pair interaction is determined 
by the environment alone. Inserting the boundary 
field n _(t) describing two opposite kinks of shape 
±2arctan{sinh[(7rc s /i)(r — Tfe)]/cos(7TXfe/L)}, k = 1,2, 
into the action @ we obtain the contribution 



G e (xi,x 2 ,T) w Kin 



cos 2 ^(zi +x 2 ) + (f^r 



(6) 



COS -f X\ COS j^X2 

which diverges logarithmically at long time scales, but 
with a weight determined by the dimensionless conduc- 



tance K of the parallel shunt. In the end, the largest 
contribution to the partition function arises from vor- 
tices nucleated at the boundary with an interaction G « 
G c = K\u(t 2 /k 2 ) and the problem maps to a system 
of charged particles in ID with fugacity A. The corre- 
sponding RG equations take the form diK = and 
d[X = (1 — K)X; in contrast to the infinite wire, the pref- 
actor K of the logarithm is invariant under the RG flow 
and we obtain a transition at K = 1 20], with a super- 
conducting phase for K > 1, while for K < 1 the vortices 
nucleated at the boundary drive the system insulating. 

In the analysis above we have determined the shape of 
the vortex-anti-vortex pair assuming Neumann bound- 
ary conditions, i.e., no currents can leave the wire. Insert- 
ing this approximate solution back into then has pro- 
vided us with the pair-interaction G ss G c = 2KIti(t/k). 
In order to check the consistency of this approximation 
we determine the correction SS due to the finite current 
if |cl>|0_ flowing through the shunt (cf. |i|). 



SS 

n 



K 
An 



du) 



\w\ |0 N 



2?r 1 + {fi/K) |coth(wL/2c s )| ' 



(7) 



with N _ w (Att/oj) sin(wr/2); performing the integral in 
shows that the corrections are indeed small at low 
frequencies. For a setup with a poorly conducting shunt 
(K <C fi) the corrections remain small at higher energies 
when t < L/ttc Si SS ~ — ^(if//z) 2 ]h(t/k). 

However, in the opposite case K 3> [i where a highly 
conducting shunt protects the superconductor the correc- 
tions are large: in the intermediate regime L/irc s < r < 
t b = KL/tt[ic s , we find SS ~ — 2K ln(r//c) + 2ix fic s T/L. 
The first term cancels the logarithmic interaction with 
the environment, while the second term describes linearly 
confined pairs. It is then appropriate to change strategy: 
for K 3> /i and high frequencies the boundary condition 
Q reduces to 0_ = and [9^0]+ = for I = and the 
phase-field for the pair is given by the L-periodic solution 
producing an interaction between the vortices 



G p (x, t) = fxln 



—) fsinh 2 
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describing defect pairs linearly confined along the T-axis 
for distances r > L/irc s , cf. Fig. [3 At the same time, 
no voltage appears over the shunt and the contribution 
from the environment vanishes, thus G ~ G p . Hence, for 
a highly conducting shunt K 3> [i the interaction starts 
with a logarithm G(r < L/irc s ) ~ ^ln[(r 2 +a; 2 )/K 2 ], goes 
over into confinement G(L/ttc s <t< tb) ~ 2tt^c s t/L, 
and ends with the low- frequency behavior G(tb < r) ~ 
K \h(t 2 / k 2 ) determined by the environment. 

In the end, we find that for all values of K the low- 
energy physics of the system is determined by the en- 
vironment. On the other hand, the renormalization 
down to low energies depends on the ratio n/K: for 
K <C /i the interaction is entirely determined by the 
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environment, while for K 3> (J. the phase-slip pairs go 
through an intermediate regime of linear confinement in- 
ducing an exponential drop of the phase-slip fugacity, 
A(t b ) ~ X(L/nc s ) exp[-G p (r B )] = X(L/nc s ) exp(-/T); 
this drop in A manifests itself in the ^-V^-characteristic. 

The above analysis results in a T = phase diagram 
as sketched in Fig. ^ The transition at K = 1 separates 
a superfluid phase for K > 1 from an insulating one at 
K < 1. In addition, we distinguish between two differ- 
ent superconducting regimes at K > 1: for /j < 2 the 
fugacity for the nucleation of small (r < L/irc s ) vortex- 
anti- vortex pairs is strongly increased as compared to the 
regime fi > 2, a leftover from the SI transition in the in- 
finite wire which is relevant at temperatures T > hnc s /L 
and drives / > 2e/j,c s /L probing these small scales. 

With the interaction between defect pairs turning log- 
arithmic on large scales t > r K = max(K, r B ), we find 
that the low-energy physics of the wire reduces to that 
of a Josephson junction with a parallel shunt R; the par- 
tition function is equivalent to that of a particle in a pe- 
riodic potential with damping rj — K/2ir |2fJ. The I— V- 
characteristic of such Josephson junctions has been stud- 
ied in detail [2(3, EJ; below, we review the main results. 
The restricted validity of the effective action S = S w +S C 
to low energy scales limits this analysis to low tempera- 
tures T < H/t k and low currents / < 2eK/irr K . 

Within the superconducting phase K > 1, the I — V 
characteristic at low drives / is calculated perturbatively 
in the vortex fugacity A (with A(r K ) — > A the renormal- 
ized fugacity at t k ). At finite temperatures the two en- 
ergy scales T and hirl /2eK define two regimes in the re- 
sponse, V ps (irI/2eK < T/h) ~ X 2 R Q I[TT K /h] 2K - 2 and 
V ps (T/H < nI/2eK) ~ (H^/eT^llT^/e}^- 1 ; the alge- 
braic characteristic at T = turns into a linear response 
at finite T and small drives I. We conclude, that in ID 
wires superconductivity survives only at T = and under 
the condition of a good protection by a high conductance 
shunt with K > 1. In addition, if K 3> /i, the linear 
confinement in the regime c s /L < r < r B produces a 
sharp voltage step Vdec ~ (^A 2 / en){L/ c s k)[ttkc s / X] 2/i_2 
at the deconfinement current /dec = 2e[ic s /L. At high 
drives I > /dec (and at temperatures T > hirc s /L) 
we start probing the wire (i.e., intervortex distances 
r < ttCs/L) and the I-V w characteristic is given by the 
result for the infinite system V pa (Trc s /L < 7r//2^e <§; 
T/h) - {\ 2 R Q I){L/c s K)[TK/h] 2 ^ z and V ps (nc s /L < 
T/h < irl/2fie) - (h\ 2 / 'en){L / 'c s k)[1 k/ e/i] 2 ^ 2 . The 
wire's I w — ^-characteristic derives from solution of the 
implicit equation V w (7 W ) = V ps (I = I w + V W /R P ). 

In the insulating phase K < 1 the phase slips de- 
scribe tunneling of the phase difference 0_ to neighbor- 
ing states <fi- ± 2-7T. In the limit K — > this tunneling 
between periodic minima is coherent and leads to the 
formation of a Bloch band e(q) of width Wo — HX/k, 
with q the quasi- momentum associated to (/>_. Apply- 
ing a small current /, a voltage 2eV w — Ti(<f>-) — d q e(q) 



is set up across the wire and all the current flows over 
the parallel shunt [2l|; the system exhibits a linear re- 
sponse V w — RI. This behavior remains valid for finite 
K < 1, but with a renormalized band width [lfil EH E3 | 
W ~ Wo[Wqt k /K\ k ' defining a linear response 

regime for small voltages 2eVw < 2eV c = max(9 9 e) ~ W. 
At higher drives, the voltage V w decreases as the phase 
slips no longer block the wire and the current flows across 
both channels. 

The SI quantum phase transition in the finite wire is 
different from the one in the infinite system; it is of the 
type inherent to finite systems coupled to a dissipative 
environment, e.g., the resistively shunted Josephson junc- 
tion or the dissipative two-state system j2^. The 
environment determines the wire's thermodynamic state 
and hence the overall shape of the Tw-Kv-characteristic, 
see Fig. ^ the wire's parameters fi and A make their 
appearance (through the voltage jump at /dec and the 
wire-dominated characteristic at large currents / > /dec) 
only in the superconducting phase K > 1 and for the 
case K 3> \i where a highly conducting shunt is pro- 
tecting the wire. The insulating phase is established at 
small values K < 1; phase-slips drive the current over the 
parallel shunt, resulting in an ohmic characteristic, un- 
til leakage through the wire sets in at V c . Here, we have 
assumed a small vortex fugacity A and hence the observa- 
tion (and consistent interpretation) of the characteristic 
voltage-peak requires a sensitive probe as V c oc A 1 ^ 1- ^. 
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